Tunneling Current Spectra of a Metal Core/Semiconductor Shell Quantum Dot 

Molecule 
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The transport properties of a metal core/semiconductor shell quantum dot molecule (QDM) 
embedded in a matrix connected to metallic electrodes are theoretically studied in the framework 
of Keldysh Green function technique. The effects of the electron plasmon interactions (EPIs) on 
the tunneling current spectra of QDM are examined. The energy levels of the QDs, intradot and 
interdot Coulomb interactions, electron interdot hopping strengths, and tunneling rates of QDs are 
renormalized by the EPIs. The differential conductance spectra show peaks arising from the plasmon 
assisted tunneling process, intradot and interdot Coulomb interactions, and coherent tunneling 
between the QDs. 



I. INTRODUCTION 

Although many of the interesting transport phe- 
nomena related to quantum dots (QDs) including the 
Coulomb blockade, 1 ) Kondo effect, 2 ) Fano resonance, 3 ) 
Pauli spin blockade, 4 ) photon (phonon) assisted tun- 
neling process 5 ) and negative differential conductance, 6 ) 
have already been extensively studied, there has been 
little reported on the effects of electron plasmon interac- 
tions (EPIs) on the tunneling current of QDs. Recently, 
EPIs have received considerable attention for their appli- 
cations in nanophotonics, biology, and the harvesting of 
solar energy. 7 ^ Much effort has been focused on the effects 
of EPIs on the optical properties of nanostructures. 8-11 ) 
Experiments can observe the strong effect of EPIs on the 
exciton spectrum of individual semiconductor QDs aris- 
ing from adjunctive metallic nanostructures 12 ' 13 ). 

In the application of QD-based biosensor, the metal 
core/semiconductor shell QDs play an crucial role. 14 ) 
Because this type QDs can be readily coupled to de- 
tected proteins through electrostatic interactions. De- 
spite the measurement of optical spectra can resolve 
molecules, 7-14 ) tunneling current spectra provide a high 
efficient means to identify molecules and address elec- 
trically a single nano-object. 15 ' 16 ) This inspires us to 
study the tunneling current through individual metal 
core/semiconductor shell QDs for the application of 
nanoscale biosensor. However, it is difficult to avoid the 
proximity effect between such type QDs since QDs are 
randomly distributed. 17 ) Therefore, we consider a metal 
core/semiconductor shell QD molecule (QDM) embed- 
ded in a matrix connected to metallic electrodes shown 
in Fig. 1 to clarify how the EPIs to influence the proxim- 
ity effect resulting from interdot electron Coulomb inter- 
actions, electron hopping process, and plasmon hopping 
between two QDs in the absence of detected proteins. 



II. FORMALISM 

Here, we consider nanoscale semiconductor QDs, the 
energy levels separation of each QD is much larger than 
their on-site Coulomb interactions and thermal energies. 



One energy level for each quantum dot is considered in 
this study. The two-level Anderson model including EPIs 
is employed to simulate the QDM junction system as 
shown in Fig. 1. The Hamiltonian of the QDM junction 
is given by H = H + H QDM + H T 



H = ^^kal^a k ^^^e k b ] kcT b k)<7 (1) 
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where the first two terms describe the free electron gas at 
the left and right metallic electrodes. a\ a (b\ a ) creates 
an electron of momentum k and spin a with energy e k 
at the left (right) metallic electrode. V k ,i,p (£ = 1,2) 
describes the coupling between the metallic electrodes 
and the QDs. d\ (di, a ) creates (destroys) an electron 
in the ^-th dot: 



Hq D m = y^ Emi, a + ^ Uini^nt^ (2) 
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where E% is the spin-independent QD energy level, and 
ne,o- — d\ a di^ a . Notations Ue and Uej describe the in- 
tradot and interdot Coulomb interactions, respectively. 
tij describes the electron interdot hopping. The Hamil- 
tonian of QD molecule described by Eqs. (1) and (2) 
has already been extensively considered for studying the 
transport properties of nanostructures. 18-20 ) The Hamil- 
tonian of electron plasmon interactions (EPIs) arising 
from the metallic nanostructure of each QDs can be writ- 
ten as Ht'. 

H t = ^ ^cJQ+y^^n^o-^+Q) — y^a;^jc]cj, (3) 
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where uj£ is the plasmon frequency of the metallic nanos- 
tructures shown in Fig. 1, and Q,£ is the coupling strength 
of EPIs. The last term involving oj£j describes the 
plasmon hopping between two metallic nanostructures. 



When metal nanostructures are close enough, the plas- 
monic modes on one nanostructure can couple with that 
on other nanostructures so that hopping of plasmons 
from one metallic nanostructure to other nanostructure 
becomes possible. 

A canonical transformation can be used to remove the 



on-site EPIs from Eq. (3), that is H n 
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where S = — J^ Vt £ n £)CF {c\ — q). 21 ) In the new Hamil- 
tonian, we have the following effective physical parame- 
ters: V£ lma = V kXa e x ^-^\ V k e 2a = V k , 2 , a e x ^- C *\ 



which is assumed to be energy- and bias-independent. 
The retarded Green function of the QD density of states 
[—IvnG r ia ) has the following expression 

oo 
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where L n is given by 
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^ je -[*i(cI-ci)-A 2 (4-c 2 )] 5 and Xi = fy/u^ i n ad- 
dition, H new has an extra term J2i^j a ^^i^iA^j + c j) 
resulting from the plasmon hopping between the metal 
nanostructures, which vanishes under free plasmon aver- 
age introduced later. Under canonical transformation, we 
note that the narrowing of QD energy level broadening 
appears and QD energy levels, electron Coulomb inter- 
actions, and electron inter dot hopping strengths become 
renormalized as well. The proximity effect between QDs 
resulting from Uf- and t\ ■ is influenced by the EPIs. 
It is possible to observe an effective negative interdot 
Coulomb interactions if we have strong EPIs and a large 
hopping matrix element of plasmon (we,j)- For simplic- 
ity, this study is restricted in the case of C/| ■ > and 
identical QDs (uo £ = uoq and \ £ = g ). 

To decuple the EPIs of H new , we take the mean- field 
average to remove the plasmon field arising from c^ — c, 
that is < e g ^~^ >= e -h 2 coth 2 {hu Q /(2k B T)) _ Baged 
on such a mean- field average, we see a reduction of 
V kAa e-^ 2 coth 2 (hu J o/(2k B T)) = v kAa X £ and interdot hop- 
ping strength t c e-^ Xl - x ^ 2coth ^ huJo ^ 2kBT ^ = t c X £ij = 
t c . This approximation is valid when the tunneling rate 
arising from the coupling between the QDs and the elec- 
trodes is smaller than the EPIs, 22 ) which is our condi- 
tion of interest. Based on such an approximation, the 
co>o-dependent tunelling rates are neglected in this study. 
Meanwhile, we assume that there is no voltage difference 
between two dots. This implies that the tunneling cur- 
rents directly involving t c can be ignored. 

Using the Keldysh-Green's function technique, 23 ) the 
tunneling currents of the QDMs in the Coulomb blockade 
regime are given by 

(4) 
where f L (R)(t) = l/[e^ e -^{R))/k B T + 1] denotes the 
Fermi distribution function of the left (right) electrode. 
The left (right) chemical potential is given by ^l(^r)- 
Pl — PR — eV a , where V a denotes the applied bias. 
Notation T denotes the equilibrium temperature of the 
left and right electrodes, e and h denote the electron 
charge and Planck's constant, respectively. Notation 
^£,L(R) = ^k V k,£,L(R)^( e ~ e *0 denotes the tunnelling 
rate from the left (right) electrode to dot 1 and dot 2, 



with a boson distribution function of Nb = 
l/( e tooo/(k B T) _ -^ and a Begsel f unction f j n ( x y 

The expressions of f < (e + nujo) and f > (e — nujo) are 
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where the effective tunneling rate of T 
F £ ,l{r)X 2 with the reduction factor of X 
exp[—^g 2 coth 2 (huJo/(2kBT))]. The dressed electron re- 
tarded Green function under H new can be obtained fol- 
lowing the procedure introduced in our previous work. 20 ) 
We have the expression of Qf a 



GiA*) = E &W( £ ) 
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where m = e - Ef + iTf/2. Y\ = (T £jL + T^ R )X 2 . The 
sum of m in Eq. (9) has over 8 possible configurations. 
We consider an electron (of spin a) entering level £, which 
can be either occupied (with probability N £ ^) or empty 
(with probability 1 — N £ ^). For each case, the level j 
can be empty (with probability a,j = 1 



Nj,a 



N, 



3,° 



Cj), singly occupied in a spin a state (with probability 
bj y a = Nj )d - — Cj) or spin a state (with probability bj, a = 
Nj,a — Cj)-, or a double-occupied state (with probability 
Cj). Thus, the probability factors associated with the 
8 configurations appearing in Eq. (9) become pi = (1 — 
Ni i& )aj,p2 = (1-N ii& )bj i& ,ps = (l-Ni^)b jj(T ,P4 = (1- 
Ni,v)cj, P5 = N ii& aj, Pq = N£ t abj t a, Pi = Ni^bj^, and 
Ps = N £j pCj. S m in the denominator of Eq. (9) denotes 
the self-energy correction due to Coulomb interactions 
and coupling with level j in configuration m. We have 
Ex = t 2 J m, E 2 = Ulj + tl/ijij - Uf), S 3 = Ulj + 
tl/iH ~ Uf e ), S 4 = 217/, i + t 2 J(N ~Uf- Uf e ), S 5 = 
U[ + t 2 J{ H - U&, E 6 = Uf + Uf d + tl/fa -Uf- Uf e ), 
S 7 = Uf + Uf d + t 2 J(nj - 2Uf e ), and S 8 = Uf + 2U e tj + 
t 2 /(/j,j — Uf — 2Uf e ). Note that the fractional numbers of 



the thermally averaged one-particle occupation number 
Ni a- and two-particle correlation functions q appear only 



on the probability weights of Q\ 



,m,a' 



Once t c = and 



g = 0, the expression of Eq. (9) is the same as Eq. (3) 
in Ref. 6. iVg j<T and q can be obtained by solving the 
on-site lesser Green's functions, 20 ) their expressions are 



N i,* = - I — 2^ F e 1 F e /m ^,m,a( e )> (10) 
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Note that £ ^ j in Eqs. (9)- (11), which are valid in the 
condition of t c /Ui <C 1. In addition, the limitation of 
< Ni ia (ci) < 1 should be satisfied. 



III. RESULTS AND DISCUSSION 

Bulk metals have very large plasmon frequencies, 
which can be ten times larger than the on-site Coulomb 
interactions of QDs. Therefore, the effects of EPIs 
on the tunneling current were ignored in the previous 
studies. 18-20 ^ The plasmon frequency of metallic nanos- 
tructures is the same order of magnitude as the on-site 
Coulomb interactions and the EPIs is strong. 13 ) One can 
expect to observe the plasmon assisted tunneling pro- 
cesses in the tunneling current spectra. To reveal the ef- 
fects of EPIs on the tunneling current spectra of QDMs, 
we initially consider the case without EPIs. The tun- 
neling current and differential conductance are plotted 
in Fig. 2 with U £ = U = 60r , Ue tj = Uj = 20r , and 
Ei = Eo-\-7]£eV a . There is a large voltage across the junc- 
tion, therefore the shift of energy level Eg arising from 
the applied bias is considered by rj£eV a . We have adopted 
r]£ = 0.5 based on the assumption that QDs are located 
in the central position between two electrodes. The black 
lines show the typical staircase structure of the tunnel- 
ing current and an oscillatory differential conductance 
with respect to the applied bias arising from the intradot 
Coulomb interactions. These structures will be washed 
out with increasing temperature. We will focus on the 
transport behavior throughout at the low temperature 
of ksT = lTo- In the presence of interdot Coulomb in- 
teractions (red lines), new staircase structures appear in 
the tunneling current. Five peaks in differential conduc- 
tance labeled from V\ to V$ result from electrons of the 
left electrode through the resonant channels of e\ = Eo, 
e 2 = E + Ui, e 3 =^ + Uo, e 4 = E + Uo + Ui, and 
€$ = Eo + Uo + 2Ui. In the presence of t c (see the blue 
lines), each peaks (Vi, V2, V4 and V5) split into the bond- 
ing (BD) and antiboding (ABD) states. Such structures 
can be depicted by using a single molecule with Eq — t c 



and Eo+t c states filled with one, two, three and four elec- 
trons. The electron filling of such a QDM satisfies Hund's 
rule. In addition, peaks V2 and V4 have extra peaks, 
which correspond to the spin singlet states (two electrons 
and three electrons). For instance the two electron singlet 
state has a resonant pole e = Eq + Uj — t 2 /(Uo — Ui) re- 
sulting from the p<i in Eq. (9), which is different from the 
two electron triplet state with pole e = Eo + Ui±t c from 
the ps of Eq. (9). The differential conductance structure 
resulting from three electrons can also be analyzed from 
the pe and pr of Eq. (9). On the basis of results in Fig. 
2, we find that the interdot Coulomb interactions play an 
important role in distinguishing between the configura- 
tions of one electron, two electrons, three electrons and 
four electrons. Many theoretical works have been devoted 
to investigate the tunneling current through parallel QDs 
for the applications of quantum computing. 24 ) Neverthe- 
less, there still lacks a comprehensive theory to reveal the 
spin states of parallel QDs. The retarded Green function 
of Eq. (9) provides a closed form expression to distin- 
guish eight configurations in the parallel QDs. 

Figure 3 shows the tunneling current (J) and differ- 
ential conductance (dG) for different strengths of EPIs 
at Uo = 60r , Ui = 20r , t c = 0, u id = oo = 20r . 
With increasing the strength of EPIs, the energy levels 
of the QDs, intradot Coulomb interactions, and interdot 
Coulomb interactions are renormalized by EPIs such as 
Ef = E £ -g 2 uJo, Uf = Ui-2g 2 uJo, and U^ = Uij-2g 2 uj . 
The current spectra and differential conductance change 
considerably. Five peaks of dG labeled from V\ to V5 at 
g=0 are shifted to the low bias regime. The magnitude 
and width of these peaks become smaller and narrower 
with increasing g. This is attributed to the current re- 
duction factor of e~ 9 ( 1+A/£? ) [see Eq. (6)] and reduction 
of the tunneling rates Y\ = TgX 2 . In addition to these 
five peaks, there is a satellite peak labeled n=l, which 
arises from one plasmon assisted tunneling process. For 
the blue line, we see two one-plasmon assisted tunneling 
peaks corresponding to e = Eq+ujq, and e = Eq + Uj+uo. 
As a consequence of the very low temperatures, the struc- 
tures arising from e = Eq — c^o, and e = Eq + Uf — cjo, 
which are contributed from the first term of Eq. (5), 
correspond to electrons of the QDM with energy levels 
Eq and Eq + Uf to have one plasmon emission process 
to escape out the QDM. This processes are obviously 
suppressed in the small bias regime resulting from small 
electron population of the QDM. 

Since we consider t c = 0, the spin degree of freedom 
can not be resolved in Fig. 3. 6 ^ To further understand 
the EPIs on the current spectra of QDMs with finite £ c , 
we plot the tunneling current (J) and differential con- 
ductance (dG) as a function of the applied bias for the 
case of Ui = 100r , U £J = 40r , and t c = 6T Q in Fig. 
4. The first two peaks indicated by the black line for 
dG correspond to the BD and ABD states. The follow- 
ing three peaks are similar to those indicated by the blue 
line in Fig. 2. They result from the spin singlet and triplet 
states. For g = 0.5, these spin-dependent spectra of dG 



just shift to the low bias regimes, whereas the separa- 
tion between the BD and ABD peaks in the triplet state 
is not changed. The exchange energy of the singlet state 
—t^/{Uf — Uf •) is also invariant. For g = 1, the tunneling 
currents are enhanced resulting from the absence of the 
interdot Coulomb blockade. The spectra from the two 
particle states merge into that of one particle at g = 1. 
Consequently, the spin-dependent spectra of dG is sup- 
pressed. In addition, we observe multiple plasmon as- 
sisted tunneling processes (n=2). 



strong EPIs coupling. Because of the hopping of plas- 
mons between two metal nanostructures, the indirect in- 
terdot plasmon mediated electron electron Coulomb in- 
teractions appear. We predict that the multiple plasmon 
assisted tunneling processes can be observed in the tun- 
neling current spectra of metal core/shell semiconductor 
QDs for strong EPIs. In the presence of detected pro- 
teins, which will glue to the QDMs, the DOS of QDMs 
is changed. Therefore, the measured tunneling current 
spectra are tilted to judge the identity of detected pro- 
teins. 



IV. SUMMARY AND CONCLUSIONS 

In this study we analyzed the effects of homogenous 
EPIs on the tunneling current spectra of QDMs in the 
absence of detected proteins. As a result of the renormal- 
ization of the energy levels of QDs, intradot and interdot 
Coulomb interactions, and the tunneling rates, there is 
a significant change in the tunneling current spectra for 
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FIG. 1: A metal core/semiconductor shell quantum dot 
molecule embedded in a matrix connected to metallic elec- 
trodes. 
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FIG. 2: (a) Tunneling current J and (b) differential con- 
ductance dG — dJ/dVa as a function of the applied bias 
for Ut = Uo = 60r , Et = E F + 30r , k B T = ir , and 
Yg^L — Yi^R — 0.5ro in the absence of electron plasmon in- 
teractions. Black lines (Uij = 0), red lines (Uej = 20Yo and 

tc 



0), and blue lines (Uej = 20Yo and t c 
the tunneling current J is in units of Jo 



= 6r ). Note that 
eYo/h. 
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FIG. 3: (a) Tunneling current J and (b) differential conduc- 
tance dG as a function of the applied bias in the presence of 
electron plasmon interactions. Black lines (g = 0.0), red lines 
(g = 0.25), and blue lines (g = 0.5). We have a plasmon 
frequency cjo = 20Fo. The other physical parameters are the 
same as those for the red lines in Fig. 2. 
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FIG. 4: (a) Tunneling current J, and (b) differential con- 
ductance dG as a function of the applied bias for different 
g values. Black lines (g — 0), red lines (g = 0.5), and 
blue lines (g = 1). We have adopted the following physical 
parameters:^ = 100r , U £J = 40r , t c = 6r , k B T = ir , 
and Fi, L = Tz, R = 0.5r . 



